Distribution of lipids in non-lamellar phases of their mixtures 

Xiao-jun Li and M. Schick 

Department of Physics, Box 351560 

University of Washington, Seattle 98195-1560 

(February 7, 2008) 

We consider a model of lipids in which a head group, characterized by its volume, is attached 

^^ . to two flexible tails of equal length. The phase diagram of the anhydrous lipid is obtained within 

^^ ' self-consistent field theory, and displays, as a function of lipid architecture, a progression of phases: 

body-centered cubic, hexagonal, gyroid, and lamellar. We then examine mixtures of an inverted 

hexagonal forming lipid and a lamellar forming lipid. As the volume fractions of the two lipids vary, 

^H ' we find that inverted hexagonal, gyroid, or lamellar phases are formed. We demonstrate that the 

^ , non-lamellar forming lipid is found preferentially at locations which are difficult for the lipid tails 

to reach. Variations in the volume fraction of each type of lipid tail are on the order of one to ten 

per cent within regions dominated by the tails. We also show that the variation in volume fraction 

is correlated qualitatively with the variation in mean curvature of the head-tail interface. 

o 

C^ ; I. INTRODUCTION 



o 



o 



0^ 



I 



o 



X 



The lipid bilayer which provides the basic structure of biological membranes is composed of a large number of 

different lipids, of which many, on their own, form non-lamellar phases. Just what xcile these non-lamellar forming 

lipids play in the properties of membranes has been the subject of much speculationa^p. That they serve important 

^ functions is indicated by the fact that cells with common lipids regulate their composition to maintain a proper balance 

'~^, between those that form lamellar phases and those that do notu. At least two major roles for non-lamellar forming 

_. -. ' lipids have been proposed. One is that, because such lipids are characterized by tails which tend to^splay outward, 

►^ ^ their presence alters the pressure profile within a bilayer permitting embedded proteins to functionH— The second, 

^>^ • noting that the lipid balance referred to above is close to a lamellar non-lamellar phase transitionQD, posits that 

f^ ' these lipids with their splaying tails may serve to facilitate the formation of structures not unlike the non-lamellar 

CO phases, structures which are characterized by volumes which are difficult for other lipid tails to fill, or regions of 

^^ non-zero curvature. In this way they could stabilize transient fusioa, intermediates. For example, some scenarios of 

membrane fusion involve the formation of a thin cylindrical stalkLTcl created from the lipids of the opposing leaves 

of the bilayers_ Such a structure has a non-zero curvature, and creates regions surrounding it which are difficult for 

the tails to filB. By concentrating irusuch regions, non-lamellar forming lipids could make such structures, and the 

jrt ■ processes they bring about, less costlylj'E2l. In a series of experiments, Gruner and co-workers have shown that inverted 

hexagonal, Hjj, phases which are characterized by periodic regions which are difficult to fill can be stabilized by the 

addition of alkanesM. Presumably the alkanes are found preferentially in these regions and lower the free energy to 

''O form these structures. Theoretical confirmation of this idea in the somessjhat analogous system of block copolymer 

Ch , and homopolymer mixtures was provided by Matsenllj. Tate and Grunertj further showed that the addition of small 

O ■ amounts of long chain (two tails of 22 or 24 carbons) phosphatidylcholine to dioleoylphosphatidylethanolamine (two 

tails of 18 carbons) stabilized the Hjj phase relative to that produced by the addition of phosphatidylcholine with 

two shorter tails of 18 carbons. This provided further indiisct evidence for the common hypothesis that difficult to 

fill, "frustrated" volumes and/or regions of high curvatureEj should be correlated with a density difference between 

lamellar- and non lamellar-forming lipids. However there appears to be neither experimental evidence nor theoretical 

C^ ' calculations that directly bear on this hypothesis, or provide an indication of the magnitude of the variation in density 

of the different lipids. It is the purpose of this paper to demonstrate and to quantify, in a model lipid mixture, this 

variation of density of the lamellar- and non lamellar- forming lipids in inverted hexagonal, Hn, and gyroid (Ia3d) 

phases. 

Our paper is organized as follows. In section II, we present the model of the lipids and develop the self consistent 
field theory in real space. In section III we develop the self consistent field theory in Fourier space. In section IV(A) 
we present our results for the phase diagrams of the single lipid as a function of its architecture, and for the anhydrous 
mixture of a lamellar forming and a non-lamellar forming lipid. We also show that there exists an "effective single lipid 
approximation" by which the latter phase diagram can be obtained from the former. In section IV(B) we present our 
results for the variation of the densities of the two lipids in the lamellar, inverted hexagonal, and gyroid phases. This 



variation is on the order of 1 to 10%. We also compare qualitatively this variation with that of the mean curvature 
of the structures. 

II. THEORY: REAL SPACE 

The model which we employ has been presented elsewhereHa, so we will be brief here. We consider an anhydrous 
mixture of ni lipids of type 1 and n2 lipids of type 2. Below we shall choose their architecture so that type 1 
lipids form lamellar phases while type 2 lipids form Hn phases. All lipids consist of the same head group of volume 
Vh and two equal-length tails. Thus we model mixtures of lipids drawn from a homologous series, such as the 
phosphatidylethanolamines studied by Seddon et al.ES. Each tail of lipid 1 consists of A^i segments of volume Vt, 
while those of lipid 2 consist of N2 such segments. For convenience, we denote Nl = NaL for L = 1, 2. The tails are 
treated as being completely flexible, with radii of gyration Rq^l = (-/Vi,a^/6)^/^ for each tail. The statistical segment 
length is a. The configuration of the Z'th lipid of type L is described by a space curve r;_i(s) where s ranges from at 
one end of one tail, through s — ol/^ at which the head is located, to s = aL, the end of the other tail. The system 
is completely described by the dimcnsionless densities of the head groups $,j (r) and of the tail segments, $j (r), 
which can be written 



$(^)(r)=z;,^5(r-r,iK/2)), (1) 

(=1 

$|^\r) ^VhY^ r S (r - n.Lis)) ds, (2) 

where Vh has been chosen as a convenient volume to make all densities dimcnsionless. The number density of tail 
segments of type L is {2N/vh)^\ , and their volume fraction is {2Nvt/vh)^t = It^t ■ The sole explicit interaction 
in the system is between head and tail segments, and this interaction energy E takes the form 

_L^[$a) +|.(^),$a) +1.^)] . ^/ [$«(r) +4f (r)] [$r^(r) + $f^(r)] dr, (3) 

where x is the strength of the interaction, and T is the temperature. The effect of a hard-core repulsion between all 
elements of the system is accounted for apDroximately by requiring that the system be incompressible. The grand 
canonical partition function of the system isEII 

m n2 p "1 "2 

^ = E 7^ / n ^^'.1 n ^^"'^ exp(-S/fcT)(5(l>W + ^^1,(1) + $(2) + ^^1,(2) _ 1)^ (4) 

The delta function in the above enforces the constraint of incompressibility. The notation J 'Dti^l denotes a functional 
integral over the possible configurations of the Z'th lipid of type L and in which, in addition to the Boltzmann weight, 
the path is weighted by the factor P[r/.L(s); 0, q;l], with 



7'[r, si,S2] = A/'exp 



8RI A, ^'1 ds I . 



(5) 



where A/" is an unimportant normalization constant and Rg = (Aa^/6)^/^ is the radius of gyration of a tail of length 
N. 

We note that because of the choice of this weight function, V, and the lack of any explicit interaction between 
chain segments to prevent their intersection, the behavior of the chains is Gaussian. This is appropriate because 
we view the chains as forming an incompressible Wfl^. Under such conditions, a flexible, interacting, polymer chain 
behaves as an ideal, and therefore Gaussian, onellS^EJ. It is not clear, of course, that an ideal lipid chain, which is 
certainly not flexible, can be treated as Gaussian. This approximation must overestimate the entropy of the tails, 
whose fewer tiieMially accessible configurations are presumably modelled more accurately by the Restricted Isomeric 
States ModeBEil. But how serious is this overestimation, is also unclear. Ultimately the efficacy of our model in 
capturing the behavior of lipids can be judged only by a comparison with experiment. This was done for the lipid 
phase behavior in referenceE3, and the comparison was very good. Even the variation with solvent concentration 
and with temperature of the characteristic period of the lamellar and hexagonal phases agreed very well. It is this 



agreement with experiment of the calculated phase behavior which provides the support for applying the model to 
the calculation of other properties, such as the distribution of lipids in mixtures investigated here. 

To proceed, we make the partition function of Eq. H more tractable by introducing into it the identity 
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where the integration on VF^ extends up the imaginary axis. We also insert such an identity for the density $J (r), 
and a similar representation for the delta function which enforces the incompressibility constraint. The partition 
function becomes 



■^ L=l 



(7) 



where the grand potential Q is given by 
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+ £;[$« + <i>(f ) , $W + $f )] -^ j dv S(r) ($« (r) + -iMp (r) + <i>f (r) + 7,<l>f ^ (r) - 1), 
and Ql is the partition functions of a single lipid, of type L, in the external fields Wj^ and W^ ; 
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To this point, no approximations in the evaluation of the partition function of the model have been made, but it 
has been put in a form in which the self-consistent field (SCF) approximation appears naturally. The need for an 

approximation arises from the fact that the functional integrals in Eq. over Wf^ and W^ cannot be carried out. 
The SCF consists in replacing the exact free energy, — fcTlnZ by the extremum of fi. We denote the values of the 
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<i>J , Wj and S which extremize VI by lower case letters. They are obtained from the following set of 



self consistent equations: 
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, so that henceforth we shall drop the superscripts on the fields Wh and wt- It 



is convenient, further, to introduce the total headgroup and total tail densities 
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and to measure chemical potentails relative to that of lipid 1. This has the effect that zi = 1, and Z2 ~ z, the chemical 
potential of lipid 2 relative to that of 1. The free energy in this approximation, flscf is 
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with E given by Eq.pl There remains only the calculation of the partition function of the single lipid I of type L in 
the external fields Wh and Wt- One defines the end-segment distribution function 



g^^Hr,s)= fvri{s)6{r-ri{s))expi-f dt 



1 . dvijt) 
SRj ' dt 



+ Wh{ri{t))S{t-aL/2) + wt{ri{t))] }, (20) 



which satisfies the equation 



ds 



with initial condition 



= 2i?2v2gW(r,s) - K(r)(5(s - aL/2) + Wtir)W'^Hv, s), 



g(-E')(r,0) = 1. 



The partition functions of the two types of lipid are, then, 

Ql= /"drgW(r,aL). 

It then follows from Eqs. |l^ and ^ that 
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where zi = 1 and zi = z. The self-consistent equations |lO| to g^ can now be solved in real space. For the periodic 
phases in which we are interested, such as the lamellar, inverted hexagonal, gyroid, etc., it is more convenient to do 
so in Fourier space. 



III. THEORY: FOURIER SPACE 

Because the densities, fields, and end point distribution function depend only on a single coordinate, they reflect 
the space group symmetry of the ordered phase they describe. To make that symmetry manifest, we expand all 
functions of position in a complete, orthonormal, set of functions /i(r),i = 1,2,3..., which have the desired space 
group symmetryE3; e.g. 



^/'W '^^'i>h,i!i{^). 



(26) 

S,,,^^JdrMr)f,{r), (27) 

where V is the volume of the system. Furthermore we choose the /j (r) to be eigenfunctions of the Laplacian 



V'/^W = -;^/Kr), 



(28) 



where D is a, length scale for the phase. We set /ij^ 1. Expressions for the unnormalized basis functions for all 
space-group symmetries can be found in X-ray tablesS. The self-consistent equations in Fourier space become 



6i,i = (ph,i +lt4't,2- 



(29) 
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To obtain the partition functions and densities, we proceed as follows. For any function G(r), we can define a 
symmetric matrix 



{Gh^^ j h{v)G{v)f,{v)dv 



(32) 



Note that {G)ii — {G)ii = Gi, the coefficient of /i(r) in the expansion of G{r). Matrices corresponding to functions 
of G{r), such as 



{e^)^^^^j f^{r)e^i^)f^[r)dv, 



(33) 



are evaluated by making an orthogonal transformation which diagonalizes {G)ij. The densities are obtained from the 
end point distribution function. After Fourier transforming the diffusion equation, Eq. El^, one readily obtains the 
solution 
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where the elements of the matrix A are given by 
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From this, the Fourier amplitudes of the densities follow from eqs. p3 and p5|; 
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The grand potential within the self-consistent field approximation, Eq |l^, becomes 
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This free energy still depends on D, the length scale of the periodic phase, which must be determined by minimization 
of the free energy. After this is done, we compare the free energies obtained for phases of different space group 
symmetries, and thus obtain the phase diagram of the system. 

The infinite set of self-consistent equations E9^ to p^ must be truncated to be solved numerically, and we have 
employed up to 50 basis functions in our calculations for the lamellar, hexagonal, and b.c.c. phases, and 100 for the 
gyroid phases. 



IV. RESULTS 
A. Phase Diagrams 

We begin by presenting, in Fig. 1, our results for the phase diagram of a single lipid with tails each of N units 
as a function of its architecture. We have defined an effective temperature, T^ = l/2xN. Note that T^ is length- 
dependent, i.e., th£_.value of T^ differs for lipids of different length N even at the same physical temperature T. In 
our previous workli£l, we found that for the system of dioleoylphosphatidylethanolamine a T^ of 0.06 corresponded 
to a physical temperature of approximately 20°C. The architecture is characterized by the volume of the head group 
relative to that of the entire lipid, / = Vh/{vh + ^Nvt). We have examined the common lipid phases, including 
the bi-continuous double-diamond phase (Pn3m), for stability, and found in addition to the lamellar phase only the 
normal and inverted versions of the body-centered cubic (ImSm), hexagonal, and gyroid laSd, phases to be stable. 
One sees that the diagram is quite reasonable and illustrates quantitatively the qualitative ideas of structure being 
driven by packing considerationscil. 

We now turn to a mixture of two lipids. We have chosen Ni — N and N2 = 1.5N (i.e., ai = 1, and a2 — 1.5) 
so that the extended length of the tails of lipid 2 are 1.5 times those of lipid 1. Further we have taken jt = 
2NiVt/vh = 2.5. With these parameters, the volume of the head group relative to that of the entire lipid is, for lipid 
1, /i = 1/(1 -I- ai-ft) — 0.2857, while that of lipid 2 is /2 = 1/(1 + a27t) = 0.2105. For comparison, the relative head 
group volume of dioleoylphosphatidylethanolamine calculated from volumes given in the literatureEj is / = 0.254. 
From the phase diagram of Fig. 1, one sees that lipid 1 forms a lamellar phase, while lipid 2 forms an inverted 
hexagonal, Hjj, phase. 

The phase diagram of the lipid mixture is shown in the solid lines of Fig. 2 as a function of the volume fraction of 
lipid 1, & = 4>f^ + ^t4'\ and the reduced temperature T* = l/2xNi. Here A^i is a constant, the length of the tails 
of lipid 1 , so that the definition of the temperature T* is independent of the composition of the mixture as it should 
be. Small regions of body-centered cubic phase near the transition to the disordered phase have been ignored in the 
phase diagram. We are unaware of experimental phase diagrams from anhydrous mixtures of lamellar- and hexagonal- 
forming lipids with the same headgroup, as we have calculated here. However phase diagrams have been obtained for 
mixtures of lamellar-forming phosphatidylcholine and hexagonal-forming phosphatidylethanolamine. Those obtained 
for anhydrous mixtures of dilinoleoylphosphatidylethanolamine and palmitoyloleoylphosphatidylcholinacJ, and for 
mixtures of dioleoylphosphatidylcholine and dioleoylphosphatidylethanolamine and 10% water by weightEII are quite 
similar to Fig. 2 as they each show a significant region of cubic phase between the inverted hexagonal phase, which 
dominates at low concentrations of the lamellar-forming lipid, and the lamellar phase, which dominates at high 
concentration. We also observe that about 20% of one lipid added to the other is sufRcient to bring about a change 
of phase, which is in accord with experimented. Lastly, if the temperature is not too low, a decrease of the volume 
fraction of the lamellar-forming lipid causes the gyroid phase to be stable at increasing temperatures. Therefore the 
addition of longer non- lamellar lipids to the mixture stabilizes the non-lamellar phase, as in the experiment of Tate 
and GrunerEj. 

We now show that there is a simple "effective single lipid" approximation by which one can obtain rather well the 
phase diagram of the mixture, Fig. 2, utilizing only the information from the phase diagram of the single lipid, Fig. 
1. For this purpose, we must find a relationship between the coordinates (/, T^) in Fig.l and (8, T*) in Fig. 2. Note 
that whereas the temperature scale T* is independent of the composition of the mixture, the appropriate scale of Fig. 
1 would vary from T^ = T* /1.5 when the mixture contained lipid 2 only to T^ = T* at the other extreme. 

To obtain the desired relationship, we first note that from the definitions of / and 7^, it follows that 2Nvt/vh — 
(1 — /)// and 2NiVt/vh = 7t, from which one obtains T^/T* = Ni/N = fjt/{^ — /)• Second, as / is defined as the 
volume fraction of the head group for a single lipid, it is natural to assume / = 8/1 + (1 — 8)/2, the volume fraction 
of the head groups, in the mixture. Hence, in the "effective single lipid" approximation the coordinates of a point on 
the phase diagram of the mixture may be obtained from those of the single lipid according to 

= -f^, (41) 

h - h 

1 - f 

T* = i-T^r. (42) 

fit 

The results of this "single effective lipid" approximation is shown by the open diamonds in Fig. 2. The approximation 
is obviously very good. Thus we can obtain very easily the phase diagram of any anhydrous mixture of our model 
lipids from the results for a single lipid, given in Fig. 1. 



B. Distribution of Lipids 

We now turn to the central results of this paper, which is the distribution of the different lipids in the various 
phases. It is convenient at the outset to define two local order parameters, ip{r) and C(r): 



,/,(r)=7/'(i)(r) + V^^Hr), (43) 

^W(r)^cb^^\r)-^^^, L = I,2, (44) 



Each order parameter, ip^^> (r) measures the local difference in head and tail segments of type L normalized such that 
the integral of the order parameter over the unit cell vanishes. Thus ip{r) measures the local difference of all head and 
tail segments. It provides information on the separation of the lipid heads and tails. On the other hand, the order 
parameter C(r) measures the difference in local fractions of the tail segments belonging to the two different lipids. 
The brackets in its definitions denotes an average over some suitably defined tail region in which the order parameter 
ipir) takes on negative values. For clarity, we employ somewhat different definitions of the tail region for the different 
phases. The average value of C(r) is zero over the defined tail region. 

We begin with the lamellar phase and show in Fig. 3 the density profiles at a temperature T* ~ 0.0454 and volume 
fraction of lipid 1, = 0.6588, a point very close to the phase boundaries between the lamellar phase and the inverted 
gyroid, G//, phase. In Fig 3(a) we have plotted the volume fractions of the heads of lipids I and 2, (j)f^ and (j)}^ , 

and the volume fractions of the tails, 7(0^ , and jt't't ■ Although the lamellar forming lipid I dominates, one can 
easily see that the density of lipid 1 tails decreases near the center of the lamellae by the order of 5 per cent from its 
maximum value of about 0.6. As the system is incompressible, this implies that the volume fraction of the tails of the 
non-lamellar forming lipid increases by about 8 per cent over the same range. To illustrate the relative difference in 
their volume fractions, we plot in Fig. 3(b) the order parameter C(x) where the tail region is defined here as ip^x) < 0. 
Again one sees that the relative change in volume fractions is less than, but the order of 10 per cent, with the longer 
lipid 2 predominating in the center of the tail region. This comes about for two reasons. First, the tails of lipid I 
are shorter than those of lipid 2, so one expects there to be less lipid I near the center. Second, the effect of the 
temperature, which is to shorten the average end to end distance of the tails, is greater on the tails of lipid I than on 
the tails of lipid 2 because T^ > T'^ . We have not plotted the relative variation of the head group volume fractions 
because it is so small, of the order of 0.5 per cent. Recall that all the head groups arc identical. 

We now turn to the inverted hexagonal, Hjj, phase. We consider the system again at T* — 0.0454 but now at 
volume fraction of lipid 1 of = 0.27 where the Hjj and gyroid phases are almost in coexistence. In Fig. 4, we 
plot the order parameter '0(r). The dark regions correspond to positive values of the order parameter where the head 
groups dominate and the lighter regions where the tails dominate. The maximum value is 0.719 and the minimum 
value —0.286. Each gradation represents a change of 10 per cent. To make manifest the difference in densities of the 
tails from each lipid, we again look at the order parameter C(r) which is shown in Fig. 5. Here the tail region we 
have averaged over is the locus of points for which the order parameter i/'(r) < —0.282. Regions with i/'(r) > —0.282 
are simply shown in white. Again the darker regions correspond to positive values of this order parameter where the 
tails of the lamellar forming lipid are relatively more probable to be found, while the lighter regions denote negative 
values where the non-lamellar forming lipid is more probable. One see that the latter lipid is more likely to be found 
in the next-nearest neighbor direction, in the region which is most difficult for the tails to fill. The maximum value 
of ^(r) is 0.0352, and the minimum value is —0.0377 indicating that the relative volume fractions vary by the order 
of 3 per cent over the tail region. To make the variation in densities even clearer, we plot in Fig. 6(a) and 6(b) 
the volume fractions of the tails of the lipids as measured along the boundary of the Brillouin zone with the angle 
6 — corresponding to the direction of the nearest neighbors, and ^o ~ ""/S corresponding to the direction of the next 
nearest neighbors. One sees that the density of the minority lamellar forming lipid is reduced by about 2 per cent in 
the next-nearest neighbor direction, while that of the majority non-lamellar forming lipid is is increased by about 0.6 
per cent. Again, the variation with direction of the density of head groups is negligible. 

In order to compare these densities changes with the curvature of the inverted hexagonal phase, we plot in Fig. 
6(c) the mean curvature along the locus of points defined by ip(j) = Oj *-e., on the curve on which the difference in 
volume fractions of all heads and tails vanishes. The mean curvature as a function of angle is well fit by 
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H{9) =Ho + ^Hn cos{6ne), (46) 



with Ho = 1.55800, Hi = -3.33 x lO'^, H2 = 3.2 x 10"'', H3 = -5.3 x 10-'^,H4 = 3.59 x lO'^, and F5 = -8.2 x IQ-^. 
One see that the variation in the mean curvature is an order of magnitude smaller than that of the variation of the 
tails of the minority lipid, and also smaller than that of the majority lipid. There is also much more structure in 
the mean curvature. It is not surprising that this structure near the region where the head groups and tails meet is 
washed out in the region of the other end of the tails. 

We now consider the inverted gyroid phase. We show in Fig. 7 and 8 two views of this phase at T* = 0.045 and 
Q = 0.4921. The phase consists of two sublattices of tubes filled for the most part with head groups. The sublattices 
are related by mirror reflection. We have chosen to plot the surfaces defined by the order parameter '0(r) = 0.5 rather 
than ■(/'(r) = for in the latter case the tubes are much more difficult to recognize. Fig. 7 shows a view along the 
[001] direction while Fig. 8 is viewed from point (0.9,-2.4,2). These figures define the coordinate system which we use 
and show a cubic cell of side unity. Were the tubes shrunk to lines, three such lines would meet at nodes, and there 
would be 16 such nodes in the unit cell shown. 

We wish to show what the relative densities of the two lipids are at the positions which are the easiest for the tails 
to fill, and those which are the most difficult for the tails to fill. These correspond to the points midway between 
nearest neighbors and next-nearest neighbors in the Hn phase. These positions are not obvious in the gyroid phase, 
but can be calculated. An example of a longest distance is shown in Fig. 8 by a solid line. The coordinates of the 
point shown, which is furtherest from any tube, is (1/2, 3/16, 5/8) while the points on the tube centers closest to it 
are (15/32, 7/32, 7/8) and (17/32, 7/32, 3/8). Thus this longest distance is diong = (66)^/2/32. The shortest distances 
link two three-fold coordinated sites on the different sublattice, and one such shortest distance is shown in Fig. 8 
by a dotted line. Two particular points on the tube centers connected by this distance and which are shown in the 
figure have the coordinates (5/8, 3/8, 7/8) and (7/8, 1/8, 5/8). The distance between these points is twice the shortest 
distance between a point midway between tubes and the tube centers themselves, and is therefore d^hort = 3^"/8. 
Interestingly, the ratio of longest to shortest distances in the gyroid phase diong/dshort = (33/8)^^^/(3)^/^ is only 
about 2 per cent larger than the corresponding ratio in the hexagonal and b.c.c. phases which is 2/(3)^/^ in both 
cases. At the temperature and composition selected, we find that at the point nearest to neighboring tubes and 
which requires the least stretching of the lipid tails, the volume fraction of the tails of the lamellar forming lipid 1 
is 7t(/)( = 0.448 and that of the non-lamellar forming lipid 2 is 74(^4 = 0.537. At the position requiring the most 
stretching of the tails, we find the volume fraction of the tails of lipid 1 has decreased to 0.431, while that of the 
longer lipid 2 has increased to 0.561. 

In Fig. 9 the order parameter C(r), which shows the relative difference between the volume fractions of the tails of 
lipid 1 and lipid 2, is shown in a cut through the gyroid taken in the [110] direction and which passes through two 
points, (—1/8, —1/8, —1/8) and (1/8, 1/8, 1/8), which require the least stretch. Only the portion of the tail region 
defined by ■!/'(r) < —0.2 is shown with any gray scale variation. The solid line connects the tube centers which lie 
closest to one another. The point midway along this line is that most easily reached by the lipid tails. The fact that 
the region around this point is dark indicates that the shorter lipid 1 is relatively more probable here. The maximum 
(black) and minimum values (white) of C(r) are 9.26 x 10^^ and —6.64 x 10^^ indicating a variation of relative volume 
fraction in the tail region shown which is on the order of 8 per cent. Figure 10 shows a cut in the [801] direction 
which passes through the point which is furtherest from the center of any tube. Again the region shown in gray scale 
variation is the tail region defined by ip{r) < —0.2. The maximum and minimum values of (^(r) are 10.35 x 10~^ 
and —6.99 x 10~^ indicating a similar variation in volume fraction. The point at the center of the bent dark line is 
furtherest from the center of any tube. The black line ends on the centers of the two tubes closest to it. The lightest 
part of the generally dark tail region is very close to the point furtherest from any tube center and indicates that the 
relative probability of finding the tails of the non-lamellar forming lipid is large there. That the lightest part does 
not correspond exactly to the point furtherest from the tube center is probably due to the fact that this point is not 
necessarily the furtherest from the head-tail interface, and thus not necessarily the most difficult for the tails to reach, 
although it is probably quite close to being so. 

In Fig. 11 we show the surface defined by the order parameter ?/;(r) — —0.2 in the region of the three-fold connectors. 
On that surface the order parameter C(r) is plotted. The largest value, which is black and at the center of this region, 
is 3.25 X 10~^ showing that the lamellar-forming lipids are more probable there. The smallest value, shown in the 
lightest gray, is —1.82 x 10"^ indicating a few per cent variation in this region. For comparison, the mean curvature 
on the surface defined by ip{r) = is plotted in Fig. 12. The smallest value of the mean curvature, 0.20, occurs at 
the center, while the largest values, 0.52, occur away from it. As in the inverted hexagonal phase, one sees that the 
preferential location of the lamellar-forming lipids is at low curvature sites. Further there is more structure in the plot 
of the curvature on a surface near where the heads and tails meet than in the plot of the distribution of tail volume 
fractions deep within the tail region. That the region near the center oLthcse three-fold connectors is characterized 
by small mean curvature is analogous to the results obtained by MatsenO for the gyroid phase of block copolymers. 
We used two hundred basis functions for this plot. 



In sum, we have employed a lipid model which has given excellent results previousl}il3 for the phase diagram of 
dioleoylphosphatidyethanolamine in order to examine the distribution of tails in a mixture of lamellar forming and 
non-lamellar forming lipids. The two lipids are characterized by the same head group, but tails of different length. 
We have shown, both in the inverted hexagonal and gyroid phases, that the tails of the non-lamellar forming lipid 
are found preferentially in regions of the unit cell which are difficult to fill, and those of the lamellar forming lipid 
are found in the regions most easily filled. The variation in volume fraction is of the order of 1 to 10 per cent. The 
difference in lipid tail density is also correlated with the curvature of the surface which, loosely, separates the head 
groups from the tails although the former shows much less structure than the latter. It will be most interesting 
to apply this model to stalk-like structures thought to be of importance in membrane fusion to deterrpine by what 
amount the presence of non-lamellar forming lipids can lower the free energy barrier to their formationEI. This work 
was supported in part by the National Science Foundation under grant number DMR9876864. 
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Figure 1. Phase diagram of an anliydrous system of lipids as a function of the temperature T^ and the volume of the 
head group relative to that of the entire lipid, / = w„/(w„ + 2Nvt) — 1/(1 + 7*). In addition to the disordered 
(D) and lamellar La phases, there are body-centered cubic (bcc), hexagonal (iJ), and gyroid G phases. The 
subscripts / and // denote normal and inverted phases respectively. The dashed lines indicate extrapolated 
boundaries. 

Figure 2. Phase diagram of a mixture of two lipids with the same head group, but one with tails 1.5 times the length 
of the other. Lipid 1, a lamellar-forming lipid, is characterized by 1/(1 + 74) = 0.2857 while the non-lamellar 
forming lipid 2, with the longer tails, is characterized by 1/(1 + 1.57t) = 0.2105. T* = l/2xNi is a measure 
of the temperature, and O is the volume fraction of lipid 1. Solid lines result from the full self-consistent field 
calculation of the mixture, while open diamonds result from the single effective lipid approximation. Very small 
regions of b.c.c phases have been ignored. 

Figure 3. (a) Volume fractions of the head groups (/)|j and (/>|j of lipids 1 and 2 and of the tail groups 7t(/'|j and 

(2) 
Jt(f>h ii^ the lamellar phase at T* — 0.045 and volume fraction of lipid 10 = 0.6588. 

(b) The order parameter ({x) = cpl (x)/ < cpi ' [x) > —(pj. (x)/ < cpi ' {x) >, where the averages are taken over 
the region '4'{x) < 0. 

Figure 4. The order parameter V'(r) — (p^ i^) ^ n (r) + (p,^ (r) — (l/Q!2)'/'t (r) in the Hu phase at T* = 0.045 and 
O ~ 0.27. The maximum value is 0.719 and the minimum value —0.286. Each gradation represents a change of 
10 per cent. 

Figure 5. The order parameter C(r) in the hexagonal phase of Fig. 4. The averages are taken over the region 
■(/'(r) < —0.282. Regions with ?/j(r) > —0.282 are shown in white. Light parts of the tail region indicate an 
excess of non-lamellar forming lipids, while darker areas represent excess of lamellar-forming lipids. 

Figure 6. (a) Volume fraction of non-lamellar forming lipid 2 evaluated on the Brillouin zone edge from nearest- 
neighbor direction, 6 = Q, to next-nearest-neighbor direction, 6 ~ 7r/6. 

(b) Volume fraction of lamellar forming lipid 1 evaluated on the Brillouin zone edge. 

(c) Mean curvature of the surface defined by V'(r) = 0. 

Figure 7. A view of the gyroid taken along the [001] direction. The surface is defined by the value of the order 
parameter V'(r) = 0.5 which is well within the region in which the head groups dominate. 

Figure 8. A second view of the gyroid surface, defined as in Fig. 7. The longest distance between a point in the tail 
region and the axis of any tube is shown by two intersecting solid lines. The shortest distance between the axes 
of two tubes is shown by a dotted line. Such a line connects nearby regions of three-fold symmetry on the two 
different sublattices. 

Figure 9. Order parameter C(r) shown in a cut which passes through the point requiring the least stretch of the 
lipid tails. The points closest together are shown connected by a solid line. The volume fractions of the tails 
are averaged over the region defined by ip{r) < —0.2. The region in which '0(r) > —0.2 is shown in white. The 
maximum value of C(r), in black, is 9.26 x 10^^ and the minimum value is —6.64 x 10^^. 

Figure 10. Order parameter C(r) shown in a cut which passes through the point furthercst from the axis of any tube. 
The volume fractions of the tails are averaged over the region defined by ip{r) < —0.2. The maximum value of C 
is 10.35 X 10~^ and the minimum is —6.99 x 10~^. 

Figure 11. The order parameter C(r) is plotted on the surface defined by the order parameter ip{r) = —0.2 in 
the region of the three-fold connectors. The largest value, which is black and at the center of this region, is 
3.25 X 10^^ showing that the lamellar-forming lipids are more probable there. The smallest value, shown in the 
lightest gray, is —1.82 x 10~^ indicating a few per cent variation in this region. 

Figure 12. The mean curvature on the surface defined by ■(/'(r) = is shown. The smallest value of the mean 
curvature, 0.20, occurs at the center, while the largest values, 0.52, occur away from it. 
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